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Abstract. We describe 5 shock wave arising from continuous initial data 
in the case of triangular conservation law system arising from "generalized 
pressureless gas dynamics model" (e.g. |12H. We use the weak asymptotic 
method 



In this paper we investigate formation of <5-shock wave in the case of triangular 
system of conservation laws: 



where uq and vq are continuous functions defined on [02, a\\ such that vq is bounded 
and uq satisfies: 

f'(u {x)) = -Kx + b, 2: e [02,01], 
and K and b are constants determined from the continuity conditions: 
f'(U 1 ) = -Ka 2 + b, f'(U ) = -Kai + b 



(1) 

(2) 



«t + Cf(«))x = 0, 
v t + (vg(u)) x = 0, 



with continuous initial data 



(3) 




i.e. 



K = 



,0 = 




For the functions / and g we assume: 



/eC 2 ([[/ ,£/i]), geCWUM), 

/">0on [C/ ,C/i], 

.g'-/">0on [UoM 

3U G (Uq, Ui) such that g(U) = f'(U). 
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As we will see in the next section, such conditions provide appearance of admissible 
5 shock wave. 

The main result of the paper is construction of formulas which smoothly and 
globally in t 6 R represents approximate solution to the problem. We have chosen 
initial data so that such formulas correspond to the process of delta shock wave 
formation from continuous initial data. 

It is well known that if the solution to equation {!]) is discontinuous function 
then unknown function v contained in ^ can contain S distribution. Such form of 
the function v is natural from the viewpoint of applications. 

The basic difficulty lies in giving the sense to the solution containing 6 distri- 
bution. Namely, the second equation of the system contains nonlinearity generally 
implying the problem of defining the product of S distribution with Heaviside func- 
tion. 

Of course, this problem appears only if we directly substitute functions contain- 
ing mentioned singularities into the equation. In that case, the product of S and 
Heaviside function is defined using measure theory [TTl [12J, [22l [T9] . Also, one can 
use regularization of <5 and Heaviside distribution and define the product as the 
weak limit of the product of the approximations [161 H8j ■ Finally, recently it has 
been proposed [7] to define solutions of 6 type (i.e. solutions containing 6 distribu- 
tions) trough appropriate integral equalities as done in the case of L°° n L 1 weak 
solutions, thus avoiding problem of multiplication of singularities (see also [13jh 

According to all said above, we see that the problem of propagation of already 
formed S shocks has been explored rather thoroughly. 

But, the problem of S shock formation is much less studied. Even if it is studied, 
it has been always done for Riemann problem and always using vanishing viscosity 
approach. One of the first result of finding global smooth approximation to the 
problem of type (p])-(jl|) can be found in [15] for the case f(u) = u 2 /2 and g(u) = u. 
There, approximate solution u e to the first equation of the system is found in explicit 
form by using vanishing viscosity approximation and Hopf-Cole transformation. 
Then, substituting u £ in the place of u in equation the equation becomes linear 
equation in v and it is solved by the method of characteristics. 

General situation using the same approach (vanishing viscosity but with the van- 
ishing term of the form et(u,v) xx ) was considered in [14j . There, for the Riemann 
initial data author proves that system ([T]),([2]) admits approximate solutions con- 
verging to 6 type distribution. Author obtains the result by using various apriori 
estimates obtained from the equations with vanishing viscosity. 

Here, we consider more general initial conditions (continuous initial conditions), 
and also give explicit formula for approximate solution to the problem. It is very 
important when generalizing result to multidimensional situation since it provides 
us to describe analytically geometric singularities [4]. 

Other possible approach in general situation is to use Oleinik-Lax formula which 
gives the solution of the first equation [11]. Still, in this case one has to add 
additional assumptions concerning behavior of the solution in the moment of bifur- 
cation. 

Assumptions on the weak continuity of the solution are not enough [6], and 
additional assumptions on boundedness are not quite justified for 5 shock wave 
type solutions. 
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In this paper, we will use the weak asymptotic method which appeared to be 
rather effective in lots of situations involving nonlinear waves formation and inter- 
action [3 m m is] - 

We give basic definitions of the weak asymptotic method. 

Definition 1. By Ov(£ a ) C 2?'(R), a <E R, we denote the family of distributions 
depending on e £ (0, 1) and t £ R + such that for any test function r)(x) £ Cq(R), 
the estimate 

(Ow(e a ),r,(x)} = 0(e a ), e -» 0, 

holds, where the estimate on the right-hand side is understood in the usual sense 
and locally uniform in t, i.e., \0(e a )\ < C T e a for t £ [Q,T]. 

Definition 2. The family of pairs of functions (u e , v e ) = (u £ (x, t),v £ (x, t)), e > 0, 
is called a weak asymptotic solution of problem (HJ)— (J4j) if 

u 6t + (f(u E )) x = O v <(e), 
Vet + (v £ g(u £ )) x = Ov(e), 

u £ —u = 0-d'{e), v £ —v = Oj}'(e), e — > 0. 

t=0 t=0 



According to the previous definitions, an approximation constructed by the 
means of the vanishing viscosity is indeed weak asymptotic solution to equation 
(JTJ) . In the case of the quadratic nonlinearity (i.e. when f(u) = u 2 ) weak as- 
ymptotic solution is constructed in [3J, and in the case when / is arbitrary convex 
function in [S]. In both of the latter papers Cauchy problems with special initial 
data of type ([3J are considered. 

Such initial data has property that the solution of corresponding Cauchy problem 
can always be represented as linear combination of the Heaviside functions (see [9] 
and Theorem 6 below). This, in turn, allows us to use the techniques of the weak 
asymptotic method [H] . 

We stress that in [3j [9] , in the case of gradient catastrophe the solution is ap- 
proximated by the smooth function which describes interaction of nonlinear waves. 

Still, in the case of general piecewise monotone continuous initial data the prob- 
lem can be solved by the same technique as in the case of initial data (J3j. We 
describe briefly the procedure 10 . 

On the first step, we replace given (arbitrary piecewise monotone continuous) 
initial data by the new initial data which differs from the original one only in the 
neighborhood (xq — e 1 / 2 ,^ + e 1 / 2 ) of the point xq from which emanates charac- 
teristic along which the shock wave will appear the first (i.e. along which we will 
have the first gradient catastrophe). In the interval (xq — e 1 ^ 2 ,xo + e 1 / 2 ) the new 
initial data will be of the type J3J (where 02 = xq — s 1 ! 2 and a± = xq + e x ^ 2 ). The 
discrepancy caused by such replacement is of order 0{e 1 ^ 2 ) in L x {Yl) sense. 

With new initial data we can construct the solution using formulas from [9] or 
Theorem 6 below. 

In this paper, we give new formulas for the construction of the weak asymptotic 
solution to problem |T]), (J3j (which are simpler then ones appearing in [9]). Then, 
after obtaining the weak asymptotic solution to ([I]) , ([3J we substitute it in equation 
([2]) and solve problem ([2]), (j4| using the method of characteristics. We will explain 
this more closely. 
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The system we consider is triangular system. The first equation is scalar conser- 
vation law with convex nonlinearity. The problem of shock wave formation for this 
equation is solved in [Sj. It was done trough the concept of 'new characteristics'. In 
a matter of fact, generalized characteristics (see [2], Definition 10.2.1) are approx- 
imated by 'new characteristics'. But, unlike generalized characteristics which are 
nonsmooth and intersecting curves, the 'new characteristics' are globally defined 
smooth nonintersecting curves along which solution remains constant. Using this 
concept, we obtain approximating solution (in the sense of Definition [5]) to equa- 
tion ([T]). This solution is smooth and we can substitute it into equation ^ and 
differentiate according to Leibnitz rule the product vg(u) appearing there. Then, 
we can solve obtained equation using standard method of characteristics. 

According to the weak asymptotic method and the procedure previously de- 
scribed, we replace problem ([I]|4|) by the family of problems 

(5) u et + {f{u e ) x = Ow(e), 

(6) v et + (y £ g(u £ ) x = e > 0, 

{Ui, x < a 2 , 
u Q {x), a 1 <x<a 2 , +Ov{e), 
U , a 2 < x 

{Vi, X < 02, 
vo(x), ai<x<a 2 , +O v >(e). 
Vo, a 2 < x 

where Ox>' will be precised in Theorem 6. 

We expose the plan of the paper in more details. 

In Section 1 we recall necessary conditions for appearance of admissible i5 shock 
wave for system (JT]), ([2]). Then, we quote result in the framework of the weak 
asymptotic method that we shall need. 

In Section 2 we construct the weak asymptotic solution to problem dTJ) , (J3]) . 

In Section 3 we construct the weak asymptotic solution to problem J2]), (J4]). 

Finally, in Section 4 we find weak limit of the constructed weak asymptotic 
solution to problem ([IE]). 

1. Conditions for 6 shock wave appearance and some weak 
asymptotic formulas 



Consider system (fTJ, ([2]) with Riemann initial data: 
(9) u|t=o 



U h x < 0, 

U r , X>0, 



(io) 4=0 



Vu x<0, 

V r , X>0. 



Since the aim of the paper is to describe formation of 6 shock waves, we want to 
determine sufficient condition on / and g which provides 5-shock wave formation 
from initial data (J9j> , (fi"0|) . In other words, we want to determine conditions on / 
and g such that Riemann problem |T]), ([2|), ([9]), (fTQ| admits solution of the type: 
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(11) u(x,t) 



Ui, x < ct, 

U r , X > ct, 



(12) v(x ) t)=\ h x<ct > +const.-t-S(x-ct). 

I Vo , x > ct 

The solution is understood in the sense of Definition f from paper [5] . 

As the admissibility conditions for S shocks we shall use overcompresive condi- 
tions (as in [H Mini [21]): 

(13) A i (E>r,V r )<c<A i (0i,l5)> »=1,2, 
where A-j, i = 1,2, are eigenvalues of system (jSJ, i.e. 

Ai(tt,«) = /'(«), A 2 (u,u) =g{u). 
From () 1 3(1 and expressions for Ai, i = 1,2 we have: 

f'(U r ) < C < f{Ul) 

(14) g{U r ) < c < g(Ui). 
The following conditions were used in [14j : 

g'>0, f">0, f'<g. 

Still, such conditions will not necessarily give S shock even if the classical solution 
u to (fTJ), §3§ blows up after certain time. Since in this paper we are interested only 
on the S shock appearance phenomenon, we shall need more restrictive conditions. 
The conditions which we shall derive below ensures S shock wave appearance if the 
classical solution to (p}, |3]) blows up. We stress that S shock wave can arise also 
in the case of less restrictive conditions on / and g but in the special case of initial 
data. 

We proceed with deriving of necessary conditions. Initial assumption is convexity 
of the function /, i.e. /" > 0. We have to find conditions on g such that (| 14[) is 
satisfied. The following condition obviously implies (fl4| : 



(15) g(U r ) < f(U r ) < c < f'(Ui) < g(Ut). 

Since /' is increasing it is clear that it has to be U r > Ui. If we assume that 
F = g — f is increasing in the interval [U r ,Ui] and that F attains zero in that 
interval, obviously (TT5")) will be satisfied (since F changes sign on [U r , U{\). We can 
collect previous considerations in the following theorem: 

Theorem 3. Assume that the functions f,g S C 2 (R) satisfy 

(1) /" > on R 

(2) g'-f">0 on [U r ,U t ], 

(3) 3U e (U r , Ui) such that g(U) = f'(U). 

If U G [U r ,Ui] then Riemann problem (OP,®, i f-? 5)) . (0) admits S type solution of 
the form IW) . 



Next, we give very important theorem in the framework of the weak asymptotic 
method (sometimes called nonlinear superposition law): 
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Theorem 4. [9] Let9i £ (x) = uji(x/e), i = 1,2, where lim u>i(z) = 1, lim = 

2 — > + 0O Z — > — OO 

arte? £ <S(R) where 5(R) is </ie Schwartz space of rapidly decreasing func- 

tions. For the bounded functions a, 6, c depending on (x, t) G R + x R we have 

(16) /(a + Wufai-aO+cMva-aOH 

/(a) + U (<P\ ~ x) {f (a + b + c)B x + f(a + b)B 2 - f(a + c)B 1 - f(a)B 2 ) + 
M<^2 - x) (/(a + b + c)B 2 - f(a + b)B 2 + f(a + c)B 1 - f(a)B 1 ) + O v >{e), 
where for p G R we have 

(17) B 1 (p) = J u} 1 (z)cu 2 (z + p)dz and B 2 {p) — J uj 2 (z)uj 1 (z - p)dz, 
and 

B^p) + B 2 (p) = 1. 

2. Weak asymptotic solution to flU, © 

First, we determine the function u £ . In the sequel we use the following notation 
(as usual x G R, t G R + ): 

ui=ui(x,t,e), Bi=Bi(p), <pi = <pi(t,E), 

9 i£ = 8 te (Lpi - x) = uji(— — -), 

e 

x d a ( \ d ~ x \ • 1 o 

die = --T^ieVPi -X) = --r-Wi{ ), I = 1, 2, 

dx ax e 

f'(U)t + a 2 - f'{ul)t-a 1 Mt) 

e e 

a\ - a 2 



x =f(U 1 )t + a 2 = f(U )t + fll = m) - m) ■ 

The function r is so-called 'fast variable'. It is equal to difference of standard 
characteristics of equation @ emanating from a2 and aj., respectively. When we 
are in the domain of existence of classical solution to (jTJ) , Q we have t — > — oo, 
while when we are in the domain where solution to |T] ([3|) is discontinuous (i.e. in 
the form of the shock wave) we have t — ► oo. 

The point (t*,x*) is the point of blow up of the classical solution to (p}, ([3|). 

Also notice that w — lim 9i £ — 9(ipi — x) and w — lim 8i e — 5((pi — x) for Heaviside 

function 9 and Dirac distribution 5. 

First, we will describe the shock wave formation process for problem ([1]), ([3]). 
We explain the procedure we shall use before we formulate the theorem. 

It is well known problem ([T]), ^ will have classical solution up to the moment 
t = t* . The choice of our initial data is such that in the moment of blow up of the 
classical solution the shock wave will be formed and it will not change its shape for 
any t > t*. This is because all the characteristics emanating from [02, ai] intersect 
in one point (t*, x*). 

So, for t > t* we have to pass to the weak solution concept. In other words, in 
the moment t — t* we stop the time and solve Riemann problem for equation ([1]). 
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Our aim is to find global in time approximate solution to {]]), ((3J which is at 
least continuous. To do this we have to avoid intersection of characteristics. 

Natural idea is to smear the discontinuity line, i.e. to take e neighborhood of the 
discontinuity line and to dispose characteristics in that neighborhood in a way that 
they do not intersect and as e —> all of them lump together into the discontinuity 
line. Of course, this will not be the standard characteristics for problem ([1]), ([3]). 
Nevertheless, along them approximate solution to our problem will remain constant. 
Such lines we call 'new characteristics'. 

Another question that arises here is how to distribute 'new characteristics' in 
the e neighborhood of the discontinuity line. The obvious way to accomplish this 
is to distribute the 'new characteristics' uniformly in the mentioned area, i.e. in a 
way that every of them is parallel to the discontinuity line. 

Since all the characteristics emanating from the interval [02, a{\ intersect in the 
same point, roughly speaking, it is enough to find the way to dispose 'new charac- 
teristics' emanating from a 2 and a\ so that they do not intersect. 

We use Theorem ?? and 'switch' functions Bi, i = 1, 2, appearing there. 

Denote by <fi, i = 1,2, the new characteristics emanating from the points a^, 
i = 1, 2, respectively. They are given by the following Cauchy problems: 

= (MP) - Bi{p))f\Ux) + cBi(p), pi(0, e) = 01 + Ae^±^, 

(18) dt 2 
-gp%{%€) = (B 2 (p) - BtWfQUo) + cBi(p), tp 2 (0, e) = a 2 - Ae^±^, 

for large enough constant A. As we shall see later, it will be necessary to extend 
a little bit the interval [a 2 ,a{\. Therefore, we have Ae ai ~^ a2 accompanying initial 
data in (fig]). Also, in (fTg]) we define: 

(19) P = ^MzMifl , 

According to what we said above, we expect that for every t > t* it should be (since 
new characteristics should be 'close' one to another for t > t*, i.e. equal as e — > 0): 

Vl (t,e)- V2 (t,e) = 0(E), t>t*, 

and (since new characteristics should be 'close' to the discontinuity line): 

d , f(Ui) ~ /(go) n( , , ^ 

More precisely, we expect that for t > t* the expression B 2 {p) — B\{p) to be 
close to zero thus eliminating nonlinearity /' appearing in the equation of new 
characteristics (TT8l) . This means that, according to Theorem ?? we have B 2 (p) + 
B\(p) = 1, B\ and B 2 are close to 1/2, and this implies that c from (fl8|) should 
be close to 2 ^j^zff^ (Rankine-Hugoniot conditions). Actually, here we use the 
following simple observation. 

Once the shock wave is formed, it continuous to move according to Rankine- 
Hugoniot conditions and it does not change its shape along entire time axis. There- 
fore, the linear equation: 

r9n s 8u cdu J(Ui)-f(Uo) 

( 20 ) «7+o« _ = ' c = 2 FT tj ' 

Ot 2 OX U\ — Un 
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and equation ([I]) with the same initial condition 

«|t=o 



Ui, x < 0, 
U Q , x>0, 



will have the same solutions. Clearly, it is much easier to solve linear equation 
(|20p then nonlinear equation |T]). Still, the question is how to pass from nonlinear 
equation {1} to linear equation (|20[) in the domains where they give the same solu- 
tion (in the case of our initial data it will be after the shock wave formation). We 
explain briefly how we do it. 

Define 'new characteristics' as the solutions to the following Cauchy problem: 

x = (B 2 -B 1 )f'(u 1 )+cB 1 , 

(21) . ai + a 2 . , , 

mi = 0, tti(O) = u (x ), x(0) = x Q + sA(x — ), x e[a 2 ,ai]. 

Thus, ifii(t,e) — x(di,t,e), i = 1,2, where x is the solution to ([21]) . We will show 
later that it is possible to choose the constant A so that for xq £ [02, 01] every t > 
we have 

dx 

5— > 0. 

This means that 'new characteristics' indeed do not intersect which in turn means 
that there exists the solution xq of the implicit equation: 

(22) x(xo, t, e) = x. 

Bearing in mind that B\ ~ B 2 after the interaction we see that, using the new 
characteristics, we have smoothly passed from the characteristics of equation (TT]) to 
the characteristics of equation (|21[) . i.e. form equation (TT]) to equation (|2T)|). 

We formalize the previous considerations in Theorem[5j The theorem is analogue 
to the main result from [9]. The problem which we consider here, i.e. problem |T]), 
([3]) can be solved in more elegant manner (see Theorem [7] below) . Still, approach 
used in Theorem [5] can be used on the case of arbitrary piecewise monotone initial 
data. Also, Theorem [5] represents motivation for Theorem [7] 

Theorem 5. The weak asymptotic solution of problem flj), (0) has the form: 

,<pi(t, e) — x. 



(23) u s (x,t) = Uo + (u 1 (x,t,s)-U )u 1 (- 



(Ui - ui(x,t,s))u 2 ( ), 



where u?j satisfies the conditions from Theorem 5 and: 

u>i(z) — 1 for z > and Ui^iz) = for z < 0. 

The functions (pi(t,s), t G R + , i = 1,2, are given by U8\) and the function p is 
given by U9\) . 

The function u\(x, t,£) is given by 

u\{x,t,e) = u (xo(x,t,e)) 

where xo is the inverse function to the function x — x(xo,t,e), t > 0, e > 7 of 
'new characteristics' defined trough Cauchy problem &21\) . 
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Proof: We substitute anzatz |23|) into ([5]): 

(u + ( Ul (x,t,e) - U Q )u l ( Px ^ £) ~ X ) + (U 1 -u 1 (x,t,e))u 2 ( <P2( ' t ' £) ~ X )) + 
V e e J t 

(f(U + (ui(x,t,e) - f/ ) Wl (^M^) + {Ul _ Ul (x,t,e))u 2 ( Mt,£) ~ X ))) 

V £ £ J x 

= V ,{£) 

After using (|16p and Leibnitz rule, and collecting terms multiplying B\ E — 

9 2e = ^(^f^), 5 l£ = ( Wl (Jeik£h-E)^ and «5 2e - ( W2 (^*f^))^ we have: 



2rJ 



(( Ul - £/o)¥>it - B 2 (p) (/(uO - /(C/ )) - Bi(p) (/(C/i) - /(f/i + U - 5 le + 

(([/! - Ui )vj« - S 2 (p) (/(^i) - f(ui)) - B 1 (p) (f(Ut + U - ui) - /(C/ ))) <5 2e - Or- 

In the sequel we write only Bi instead of Bi(p) 1 i — 1,2. 

We rearrange this expression using the following simple formula C0\ e + D9 2e = 
(C + D)6 2e + G{9 le - 

(£ + [<*-*)/-(.,>]£)(*.-« 

((tti - Ub)vit - S 2 (f( Ul ) - /(C/ )) - Si (f(Ui) - /(C/r + C/ - tti))) + 

- «i)fPat - B 2 (/(C/r) - /(m)) - B x (/(C/i + £7 - u x ) - f(U ))) 5 2e = O v >{e). 

For an unknown constant c we add and subtract the term cB x in the coeffi- 
cient multiplying (9 Xe — 2e ) and then we rewrite the last expression in the following 
form: 
(24) 

^ + [{B 2 - B 1 )f'( Ul ) + cB 1 ]^\ (9 lE - 9 2e ) + 

Bi[^ (f(u + 4- u x ) + f( Ul ) - CUl )} (e le - e 2£ ) + 

((ui - U )(p u - B 2 (/(ui) - f(U )) - Bi (/(C/r) - /(C/r + U - Ul ))) S le + 
(([/! - ui)ip2t - B 2 (/(C/r) - /(«i)) - B 1 (/(C/r + U - m) - f(U ))) 8 2e = O v >{e). 
We put 

+ [(B 2 - B l )f'{ui) + cB x ] = 0, u 1 (x,0,e) = u o (i), ieh«i]. 
The system of characteristics for this problem reads: 

, N x = {B 2 -B 1 )f'(u 1 ) + cB 1 , 

(25) 

iii = 0, ui(0) = Uo(x ), a:(0) =x 6 [02,01]. 

The aim is to prove that characteristics defined by the previous system do not 
intersect. It appears that it is much easier to accomplish this if we perturb initial 
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(26) x=(B 2 -B 1 )(-Kx + b) + cB 1 , x{0) = x + eA [ x 



data for x in the previous system for a parameter of order e. More precisely, instead 
of (|25|) we shall consider system ([2T|) (the same is done in [9]). 

It is clear that such perturbation changes the solution of (|2"Tj) for Od'{s) since 
initial condition in (|21|) is continuous. 

We pass to the proof that the characteristics given by (|2Tj) do not intersect. From 
the second equation in (|2"Tjl it follows Ui = Uq{xq), We substitute this into the first 
equation of (|2T|) and use /'(wo(ieo)) = —Kx Q + b, x 6 [a2, ai]. We have: 

ai + a 2 ' 
2~~ 

Out of the segment [a 2 — eA ai 1 2 a2 , ai + eA ai ^ a2 ] initial function is constant and 
we define the solution u\ of problem ([21]) to be equal to XJ\ on the left-hand side 
of the characteristic emanating from a 2 — A ai 1 2 a2 and to be equal to Uo on the 
right-hand side of the characteristic emanating from a\ + A ai ~^ a2 

For the functions ipi and <p 2 as the characteristics emanating from a\ + A ai ^ a2 
and a 2 — A ai ^ a2 respectively, we have (compare to (TT5|0 

(27) Vlt = (B 2 -B 1 )(-Ka 1 +b) + cB 1 , 

(28) <f2t = (B 2 -B 1 )(-Ka 2 +b) + cB 1 . 

Now, we show how to effectively determine p given by (|19p . We apply standard 
procedure (see [SOU [9])). Subtracting from ([28]) we get: 

fa>2 - Pi)t = e ( ^ 2 ~^ - ep* - (B 2 - 

Then, passing from the "slow" variable t to the "fast" variable r we obtain (we also 
use B 2 + S a = 1): 

(29) p T = l-2B 1 (j>), P - =1. 

T r — ► — OO 

We explain the condition lim £ = 1. We have from (f7f| and ([28]) 

T — * — OO T 

p = /q 2(t7 - ugKgg - + a 2 - 01 

r 2(u — Wg)i + fl 2 — a i 

Putting t = in the previous relation we see that 

(30) £ =1. 

r t=o 

When we let e — > when t = we have r — > — oo. Therefore, from (|30|) it follows 

£ =1. 

T r — > — OO 

This relation practically means that new characteristics emanating from ai, i = 1,2, 
coincides at least in the initial moment with standard characteristics up to some 
small parameter e. Still, since r — > — oo for every t < t* (which means i?i — > 0) 
we see from (|2"T|) and (|28l) that new characteristics coincides with standard ones for 
every t < t* up to some small parameter e. 

Next, we analyze (|2"9"]l . From the standard theory of ODE we see that p — > po as 
t — ► +oo where po is constant such that Bi(p ) = B 2 (p ) = 1/2. That means that 
after the interaction, i.e. for t > t* , we have 

P = = Po + C(e) =► = ¥>2 + 0(e), £ -» 0, 
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or, after letting e — > 0, for t > t* we have shock wave concentrated at (see text in 
front of theorem for notations): 

(31) <p(t) = lira <pi(t,e) = ^(t-t*)+x*. 

Now, we can prove global solvability of Cauchy problem (f2"Tj) . 

Problem (f2Tj) is globally solvable if characteristics emanating from the interval 
[aa — Ae ai ^ a2 , a\ + Ae ai ^ a2 ] do not intersect. To prove that we will use the inverse 
function theorem. We will prove that for every t we have > which means 
that for every x = x(xo,t,e), xq G [02, a±], we have unique xq — xo{x,t,e) and we 
can write Ui(x(xq, t, e),t) = Uq(xq(x, t, e)). 

Differentiating (|26[) in xo and integrating from to t we obtain (we remind 
B 2 + B 1 = 1): 

dx /"* /"* 

(32) — = 1 + eA - K / (B 2 - BAdt 1 = l + e-K / (1 - 2Bt)dt'. 

OXo Jq Jo 

For t G [0,i*] we have (notice that 1 - Kt* = 0): 
7— = 1 + eA — K / dt + K 2B x dt > 



dxo 

1+eA- K I dt + K I 2B x dt = eA + K I 2B x dt > 0. 



So, everything is correct for t < t*. 

To see what is happening for t > t* , initially we estimate 1 — 2B\(p) when 
t — ► 00. From equation (|2U)) we have (we use Taylor expansion of B\ around the 
point p = po): 

p T = l-2B 1 (p) = -2(p- Po )B' 1 (p), 
for some p belonging to the interval with ends in p and po- From here we see: 

P - Pa = (p(to) - Po)exp( / -2B[(p)dT') = (p(-r ) - p )exp((T - t)2B[(pi)) 



for some fixed p G R and p\ G (p(to) , p(t)) C [p(T ),p ]- We remind that 
B[(pi)) > c > 0, for some constant c, since -Bi is increasing function and pi 
belongs to the compact interval [p(r ),po]. Letting r — > 00 we conclude that for 
any N G N 

p- P0 = O(1/t n ), t^oo. 
From here we have p r = C(l/r ), r — > 00, since: 

lim — — = lim 1 ~ 2B ^P) = lim -2S^(p) = -2B[(p ) = const. < 

T-*CO p — po T^OO p — po T^OO 

This, in turn, means that for every JVeN and t > t* we have 
(33) l-2B 1 (p)=p T = 0(t- n ) = 0(s n ), e->oo, 

since for fixed t > t* we have r = _> qq as £ _> 0. 
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Now we can prove that tP 2 - for t > t*. We have 



(34) 



dx 

dx 
dx, 



r 

- = 1 + sA - 2K / (1 - 2B{)dt' = 
o Jo 

1 + eA - 2K [ (1 - 2B 1 )dt' - 2K [ (1 - 2Bi)dt' = 
Jo Jt* 

eA + 4 f Btdt' - 2K f (I - 2B{)dt' > eA - 2K ( [l-2B 1 )dt'. 
Jo Jt* Jt* 

Recall that 

B 1 = B 1 (p(r))=B 1 (p(^)). 

e 

Consider the last term in expression (|34D : 

2K J\l - 2B 1 )dt! = 2K j (1 - 2B 1 (p(^p-)))dt' = 



t* <t' <t ==>• 0< z < 

^o(* 



2isTs / E (1 - 2B 1 (p(z)))dz < e2KC, 
Jo 



10 

where 

poo 

C= (1 - 2B 1 {p{z)))dz < oo, 
Jo 

since from (J33J) we know 1 — 2Bi(p(z)) = 0(^ _ -' v ), z ^> oo and ]V £ N arbitrary. 

Therefore, for A large enough (more precisely for A > C)we have > what 
we wanted to prove. 

Next step is to obtain the constant c. We multiply (|24|) by ry e Cq(R), integrate 
over R with respect to x and use (pH]) (so, we remove the first term in (|2"4l ): 

Si[^- (/(t/i + C/o - ui) + /(ui) - cui)] (6 le - 6 2e ) r ] (x)dx+ 
ax 

((«i - U )^u-B 2 - f{Uo))-B 1 C/(C7i) - /(C7i + C/ - «i))) Sie+ 

((C/i - «i)^ 2t + S 2 (/(ur) - f{U x ))+B x (/(C/ ) - /(CTi + C/ - «i))) <5 2£ = 0(e). 

We apply partial integration on the first integral in the previous expression to 
obtain: 



(35) / Bi[/(C7i + C/ - ui) + /(tti) - cui] (0 i£ - 0a e )r/(a:)da:+ 

((ui - C/o)^u ~ S 2 (/(mi) - /(C/o)) + Bi (/(mi) + /(C/o) - cm)) V (x)S ls dx+ 
((Ui - u 1 )<p 2t +B 2 (f( Ul ) - f(U 1 ))-B 1 (/K) + /(C/0 - cu l ))i 1 {x)5 2e dx = 0{e). 



Now we use the choice of regularizations u>i, i = 1,2. Since, 
(5i e = -w!(— ) = for a; < ^i 

£ £ 

1 f <p 2 -x^ 
o 2e = -u> 2 { ) = tor x > ip 2 

£ £ 



DELTA SHOCK WAVE FORMATION 



13 



we have from (13511: 



epB 1 / [f{Ui + U - ui) + f( Ul ) - cui]-^ 2 -lr,'(x)dx+ 

((«! - L/oteu - S 2 (/( ttl ) - f(U Q )) + B 1 (f( Ul ) + f(U ) - cm)) ?7(x)5 le dx+ 
^ (((7! - U i)(^ 2t + B 2 (/(m) - f(U 1 ))-B 1 (f(ui) + f{U{) - cu x )) V (x)5 2e dx = 

— oo 

(now, we use Wi = f/i for x > (pi and tii = Uq for a; < ^2) 
(36) e P B 1 [ [f{Ux + U Q - m) + /(t*i) - cit /^"^ ^^ 

B\ (2f(U )-cU )i 1 {x)8 le dx- 
Bi (-2/(£/i) + ct/i) 5 2e r)(x)dx = 0(e). 

-OO 

To continue, notice that we have \pB\\ < 00 for every r^R. Namely, 
(37) 

1 1 

\pBi(p)\ — > as t — > — 00 since in that case B x (p(t)) ~ Bi(r) ~ — ~ — tt, 
\pBi(p)\ — > pqBi(pq) as r — > oo since in that case p — > po- 
This fact, together with the fact that Si £ 5(</?j — x), reduces expression (j3"6")l to: 

(38) Bi (2/(f/ ) - cC/ ) 77(^1) - Si (-2/([/i) + cU x ) V (<p 2 ) = O(e). 
Rewrite this expression in the following manner: 

B x (2(/(C7 ) - f(Ui)) - c(U Q - E7i)) n{ipx) + B x {-2f{U x ) + cUi) (77(^2) - = 
B x (2(/(C7 ) - f{Ui)) - c(U - Ui)) V (<pi)+ 

epB l{ p) (-2/(^0 + cUl ) fetgl) =133 

Bi (2(/(t/ ) - /(EM) - c(C/ - C/i)) Jjfo) - 0(e)- 
From here, we see that the last relation is satisfied for 

(39) c=2 /(^)-/aM 

The theorem is proved. □ 

Remark 6. In the case such as our, when Uq and U\ are constants, is possible to 
replace formula ([23]) by 

U e {x,t) = u(x (x,t,e)), 

where, as before, the function Xo is the solution to implicit equation (|22[) and u are 
initial data ([3]). 

The proof of this fact obviously follows after comparing the trajectories. We 
give precise formulation in the next theorem. We leave it without proof since it is 
completely analogical to the proof of the previous theorem. 
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The difference between the previous and the next theorem is in the form of 
characteristics along which we solve our problem. 

In the previous theorem, for fixed e, the weak asymptotic solution u e to (JTJ) , 
([3]) was generator of continuous semigroup of transformations (since characteristics 
intersect along x — <pi) and in the following theorem the weak asymptotic solution 
u e to (jSJ) forms continuous group of transformation since appropriate charac- 
teristics do not intersect. Still, approach from the next theorem can be used only 
in the case of special initial data. 

Theorem 7. The weak asymptotic solution u e , e > 0, to Cauchy problem 

(40) u t + (f(u)) x = 0, tt| t =o=fi(x), 
is given by 

(41) u e (x,t) = u(x (x,t,e)), 

where xq is inverse function to the function x — x(xo,t,e), t > 0, e > 0, of 'new 
characteristics ' defined trough the Cauchy problem: 

ai + a 2 



,^ x = f'(u e )(B 2 (p)-B 1 (p)) + cB 1 (p), x(0)=x +eA [ x 
(42) V e 

u £ = 0, u £ (0) = u(xq), xq G R. 

where A is large enough, the functions B\ and B 2 are defined in Theorem ??, 
constant c is given in (20\) and p — p(ipo(t)/s) is the solution of Cauchy problem 

m 



The following corollary is obvious. It claims that the weak asymptotic solution 
defined in arbitrary of the previous theorems tends to the shock wave with the 
states Ui on the left and Uq on the right (see (f3Tjl ): 



Corollary 8. With the notations from the previous theorems, for t > t* the weak 
asymptotic solution u £ to problem {Tp, we have for every fixed t > 0: 



u e (x,t) 



Ui, x < f(t-r) +x*, 
U , x> §(*-**) +x*, 
where — v means convergence in the weak sense with respect to the real variable. 



3. The weak asymptotic solution to flU), flU) 

At the beginning of the section, we explain some general moments. 

The plan is to substitute smooth function u e given by (|2li|) into ©. Thus, we 
obtain equation ([6]). Augmented by initial data this linear partial differential 
equation of the first order has global differentiable solution. 

But, as e — > we can have discontinuities in v — lim v s not only on the line 

£—+0 

on which the shock wave of standard admissible weak solution u of (p}, Q is 
supported, and which appears for t > t* . Also, for t < t* discontinuities can appear 
along standard characteristics for problem ([T]), Q emanating from the points a 2 
and a\ . Discontinuity arises due to non-smoothness of initial data ([3]) in the points 
a 2 and ai (see Example [T3j) . 

The situation is different if instead of ([23| we put in ([2]) in the place of of u 
the function u e given by (|4ip in Theorem 7. The function u e is not smooth since 
we do not use regularizations of Heaviside function to smooth weak discontinuities 
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appearing in the initial data. Accordingly, for t < t* it admits the same behavior 
as the standard admissible weak solution u of problem fl}, (J3]) (it can produce 
discontinuities in v £ even for t < t*). 

Therefore, for the sake of consistency (we have the same situation for e > and 
as e — » 0), we will use the function given by (141 [1 in the place of u appearing in l[2|). 

Weak asymptotic solution to Cauchy problem (J4j) we will solve separately 
in five areas of (x, t) plane in which the solution is certainly smooth. Then, we will 
connect solutions in those domains and prove that the function formed in that way 
represents weak solution to our problem. 

In order to single out those domains we substitute (|4Tj) into ([2]) and, formally, 
use Leibnitz rule for derivative of product: 

(43) v £t + g(u £ )v £X = -(g(u £ ) x v £ . 

Remark 9. We say "formally" since the function v £ can have discontinuities, at 
least for e = 0. We repeat that this is caused by non-smoothness of the new 
characteristics if we use (|4Tj) as well as standard characteristics in xq = a-i and 
.t = a 2 . 

The system of characteristics corresponding to (PL3"|) . @ is: 

(44) X=g(u e ), X(0)=x , 

v £ = - v £ (g(u £ )) x , v e (0) = v(x ). 

We prove global resoluteness of this ODE system for xq G [aa — eA a2 ^ ai , a\ + 
eA a2 1 2 ai ]. According to the inverse function theorem it is enough to prove that 
along entire temporal axis we have 

dX „ 
t^>0. 

ox 

Denote by J = where x — x(xq, t, e) is the function defined by Cauchy problem 
(|26|) . We have proved in the previous theorem that J > for every t > 0. Recall 
that 

u £ (x,t) = u (x (x,t,s)), 
where iio is inverse function to the function x defined trough (|26|) . From 
have (we write below g'(uo) — g 1 (uq(xq(x, t, e)))) 



d dX . . dx dX . . dX dX 
at oxq oa oxo oxo ox 



t=0 



After integrating this differential equation with respect to the unknown function 
4P- we obtain: 

dX f* 

— = cxp( / g'iuoWoJ^dt') > 0, t > 0, 
ox Jo 

which implies existence of inverse function xo — Xa{X,t,e) along entire temporal 
axis, which, in turn, implies global resoluteness of problem (|44j) . 

Denote by tp*, i = 1, 2, solutions of the following Cauchy problems: 

X =g(u £ ), 
X{0)= ai , i = l,2. 

Now, we can introduce domains in which we will separately solve Cauchy problem 
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Remark 10. In the domains to be introduced, it is equivalent to say Cauchy problem 
(|43|) . |4]) instead of {2j), |4| since in those domains problem |2]), (j4|) has smooth 
solution. 

We set 

Di = {(x, t)\x< tp 2 }, D 2 = {{x, t)\x> ipi}, 

D 3 = {(x,t)\ ip 2 <x < ip 2 }, D 4 = {(x,t)\ ipl < x < ipx}, 

D 5 = {(x,t)\ <p* 2 <x< ifil) 

On the beginning, we prove that those domains are disjunct. Accordingly, we 
inspect relations between the functions (fi and <p* } i = 1,2. We have to prove the 
following fact for every t G R + : 

(45) ip 2 (t , e) < ip* (t , e) < y\ (t , e) < <p x (t , e) , 
First, we prove that 

(46) cp 2 < cp* 2 . 
In the moment t — we have 

(47) (cp 2 y t = f{U!){B 2 - Si) + cB 1 and {&)' t = g(U t ), 
and (see (fl5)l) 

g(Ux) > /'(f/i) > /'(C/i)(B 2 - B,) + cBx, 

since f'{U\) > c/2. Using well known theorem from ODE-s ("who goes slower does 
not reach further", [1]) from (|47|) we see that in some neighborhood of t = we have 
ip 2 < <p 2 - Assume now that to is the smallest t > such that </?2 = </?2 ■ In this case 
we have the same situation as in the moment t = 0, i.e. there exists neighborhood 
(to, to + <5) such that (p 2 < f 2 in (£oj^o + Continuing like this we see that we 
indeed have f4*6")) . 

In the completely same manner we prove that 

(48) <pl < ipx. 
It is remained to prove that: 

(49) 

This directly follows from the fact that characteristics of problem (|4"4"1) do not inter- 
sect. That means that relation between two characteristics remains the same along 
entire time axis. Therefore, 

<P* 2 = X(a 2 ,t, e) < X(ai,t, e) = tp\, 



since a 2 < a\. This proves (|491 . 

Collecting (|46|l . (|48|) and (|49|) we obtain (j45|) . Very important implication of 
relation (|43]) and the fact that for t > t* 

limine) = £(t_t*) + aj*, »=1,2 

is the following. For t > t* we have: 

(50) lim = ^(t — **) t = 1,2. 

We remind that the constants t* and x* are introduced in front of Theorem O 
Next, we solve problem (|4*5| . (|1]) separately in domains Di, i — 1, 5. 
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In domains Di and D 2 we have u e = const, and therefore the characteristics 
corresponding to v e there are straight lines. More precisely, we have: 



We use standard method of characteristics. Note that in this case characteristics 
emanate from the lines x = ifi, i = 1, 2, and not from x axis as usual The system 
of characteristics for ones emanating from the line 951 has the form: 

X =g(u e ), 

v e = - v e (g(u e )) x , 

X(t ) =<pi{t ) = x , v E (t ) = V Q . 

and for the characteristics emanating from the line ip 2 has the form: 



Global solvability of this system can be proved in the same way as one for the 
system (|33|) . 

Next step is to solve second equation of ([5"Tj) (or analogically of ([51])). We have 
for problem (|5T|) : 



v e (x,t) = Vi, (i,t)efli, 
v E (x,t) = V , [x,t)eD 2 . 



Another two domains are: 



D 3 = {{x,t)\ tp 2 <x < ^2}, D 4 = {(a;,*) I ipl < x < ipi}. 



In those domains we solve the following Cauchy problems: 



Vet + g{u e )Vex = -{g{u e ) x V e , 

Ve\x=ip x = Vo (initial data for the first Cauchy problem), 
Ve\x=ip 2 = Vi (initial data for the second Cauchy problem). 



X =g(u e ), 
i> e =- Ve(g{u £ )) x , 

X(t ) =<fi2(to) = X , V £ (t ) = Vi. 



(-g(Ue))xV e 





■t d dX 



(51) 



ax ' 

9 .TO 



and, similarly, for ([51]) : 



(52) 
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Previous implies: 



dx 

v e (x,t) = Vi— —(x,t,e), (x,t) e D 3 , 



dx, 



v e (x,t) = V :i ^(x,t,E), (x,t)eD A . 

where X03 = Xo(X,t,e) = ^i(toi) and X04 = xo(X,t,e) — <p 2 (to 2 ) (for appropriate 
toi, i = 1,2, depending on (X, t)) are inverse functions to the function X determined 
by (jnj) and (f5"Tj) . respectively. 

Finally, we solve problem ([321), © m the domain: 

D 5 = {(x,t)\ if* 2 <x< (pi). 

We apply similar procedure as in the previous case. The solution in this domain is: 



v e (x,t) = v (x (x,t,s)) :: ^-(x,t,£), 

where xos = xq (X, t, e) is inverse function to the function x determined by (|44[) (xq 
restricted on [a 2 - eA^±^, oti + eA^±^}). 
Thus, we have proved the following theorem: 

Theorem 11. The weak asymptotic solution to problem ^ is given by the 

formula: 



(53) 



V e (x,t) 



Vo, ,(x,t)eD u 
V ^(x,t,e), (x,t)€D 3 , 



v (x (x,t,e))^-(x,t,e), (x,t) e D 5 
Vi^{x,t,e), (x,t)eD 4 , 
Vx, (x,t)eD 2 . 



Let us analyze the function v £ more closely. We have to inspect its behavior on 
the boundaries of the domains -Dj, i = 1,2,3,4,5, since in the domains we know 
that equation (|43)) is satisfied. 

On the lines Di n D 3 and D 2 H D4 the function is continuous and therefore, in 
the domain Di U D 2 U D 3 U D4 the function v e represents admissible weak solution 
to PI). 

On the lines D% n -D5 = (p 2 and D 4 n -D5 = (p* the function v e can have discon- 
tinuities (see Remark [9] and Example [13]) . Therefore, we have to check if Rankinc- 
Hugoniot conditions are satisfied on tp% , i = 1, 2 for equation (|43j) given in divergent 
form: 

(54) v et + (v e g(u s )) x = 0. 

Accordingly, we have to check: 



(ft 

Since u e is continous function from here it follows: 

X — Lp* 

which is exactly the definition of the function ip* , i — 1,2. 

This concludes the proof that the function v e represents weak solution to problem 

in, ©. 
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4. Weak limit of the solution 



It remains to inspect the weak limit of the weak asymptotic solution (u e ,v e ) 
of problem @, © for t > t* (since for t < t* we have classical solution of the 
considered problem). 

We have already known from Corollary 9 that for t > t* we have: 



(55) 



w — u s (x, t) 



U 1} x < §(*-**) +x* 
U , x > §(* — **) +ac* 



So, we have to inspect weak limit of v e . ore precisely, in this section we will prove 
the following theorem: 

Theorem 12. For every fixed t > t* the function v £ given by i53\) satisfies as 
e -> 

(56) 

\V , x>§{t-t*)+x* 
Vi{a 2 + g{U x )t - £(t - f) - x*) + Vb(|(* - f) +aj* - o x - g(U )t)+ 

w (x )(ia;oJ(5(a; - ~(i - t) - x ), 

where means convergence in the weak sense with respect to the real variable. 

Proof: To begin, note that we can write function v e from (|53|) in the following 
manner: 

dx 



where 
(57) 



v s (x,t) = v(x (x,t,e)) — (x,t,e), 



'x-g(U{)t, (x,t)eDi, 

Xoi(x,t,e) ~ g(U 1 )ifi2 1 (xoi(x,t,e)), (x,t) e D 3 

(here first we go by Xq 3 to the line ip% so that Xq 3 (x, t, e) = ^2(^0) 
and then proceed to the line t = along the straight line x — g(Ui)t), 



x (x,t,e) 



= < 



' 05 



(x,t,e), (x,t) e D 5 , 



x l(x,t,s) - g{U Q )ip l 1 (x ^(x,t,e)), (x,t) G D 4 , 

(here first we go by Xq^ to the line ipi so that x^ i {x,t,e) = <£>i(io) 
and then proceed to the line t = along the straight line x — g(Uo)t) 
[x-g(U )t, (x,t)eD 2 , 



and 



dx 

— (x,t,e) = ^ 



1, (x,t)eD 1 , 
Mel)*, 

(x,t)eD A , 
1, (aj,i)GD 2 . 
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We take 77 e Cq (R) and write using (J53 



v e (x,t)r)(x)dx = / v e (x,t)r)(x)dx + / v s (x, t)rj(x)dx+ 

J — 00 J ^p2 — e 

<fil r'-Pi+t poo 

v E (x, t)rj{x)dx + / v e (x, t)rj{x)dx + / v e (x, t)rj(x)dx = 

ViT)(x)dx+ / Pi— — rj(x)dx + / v (x (x 1 t,e))— — r](x)dx+ 

-CO Jip2 — £ ® X J <P2 

~dx~ ' voV{ x )dx. 



Here, we have written <pi ± e in order to avoid possible ipi = ip*. 

Then, we use the change of variables x = X(xq, t, e) where X is inverse function 
of the function x n = xo(X,t,e) given by ([57)) . We have: 

(58) 

/IP2-E 
v E (x, t)r](x)dx+ 
-00 

£12 pai 

V 1 ri(x(xQ,t,e))dxo + / v (x )r)(X (x , t, e))dx + 

((p2—e,t,e) Ja 2 
xo(ipi+£,t,e) poo 

Vor)(x(xo,t,e))dx + / V r){x)dx, 

a-i J ipi+e 

and we remind that: 

xo((pi,t,e) = (pi + £-g(Uo)t, xo(ip 2 ,t,e) = ip 2 -e-g(Ui)t, 

and for t > t* we have (see (J27J and p8])): 

x («^i + £,*,£) -> -{t - t*) + x* - g{U )t, e -> 0, 

a; (^ 2 - £, i, e) -» |(t - t*) + x* - g(C/i)t, £ -> 0. 

Accordingly, for t > t* after letting e — > we have from (f55|) exactly (f3^|) . 
This concludes the theorem. □ 

It remains to give a comment concerning admissibility of the solution. Actually, 
it follows from assumptions on /' and g quoted in Theorem [3] providing: 

(59) g(U 1 )<f{U 1 ), f'(Uo)<g(U ). 

Admissibility of the shock wave appearing in the solution to problem |T]), (J3]) implies 

f'(Ui) < c/2 < f(U ) 
which together with (f5"9")) implies: 

g(Ui) < c/2 < g(U ), 



which proves overcompressibility of the shock and S shock wave appearing in ([55 
and (E 
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Example 13. In this example we look for the classical solution to the following 
problem: 

.1 



u t + (^w 2 ) a 
v t + (2uv) x 



-- 




for initial functions given by: 



t=o 



t=o 



1, x<-l 

—x, — 1 < x < 1, 

-1, 1<X, 

'l, |x|>l 
a; 2 / 3 , -1<x<1. 



f/smg previous notation on this case, we have: 

a\ = 1, a 2 = — 1, ?7 = 1, Ug = — 1 
b = 0, K = l. 

Also, for the inverse function xq — xq(x, t) we have: 

'x-2t, x<t-l, 
-2 t-l<x<-(t-l) 2 , 
X —(t — l) 2 < x < (t-lf, 



X<3 = < 



(i-ty 



, (t - if < x < 1 - t, 



^x + 2t, l-t<x. 

From here we can immediately compute j^j- and then 

dx 1 
dx~o~W 

The solution of the problem is: 

1, x > t - 1, 



(60) 



u(x, t) 



(61) 



v(x, t) = 



j^-, t-l<x<-t + l, 
-1, -t+l<x. 

1, X < t - 1, 

i-1 <x< -(t-1) 2 , 

(T^, -(t-l) 2 <X<(t-l) 2 , 

(t-l) 2 <x<-t+l 



k l, X>l-t. 

Notice that the function v has discontinuities on the following lines: 

x = {t-l) 2 , x = -{t-lf. 

and that on that line Rankine-Hugoniot conditions are satisfied. 

Notice, also, that if instead of x 2 / 3 in the initial data for v we would put x 2n the 
solution of our problem would be continuous for t < t* . 
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